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Введение. Рассматривается новый подход в формировании методологического базиса системного анализа в 
приложении к задачам динамики механических колебательных структур. Целью исследования является 
разработка метода оценки свойств механических колебательных систем с учетом сил вязкого трения на основе 
частотных функций и функции демпфирования, зависящих от так называемого коэффициента форм связности, 
представляющего собой отношение характеристик обобщенных координат. 

Материалы и методы. Используемые графоаналитические методы оценки динамических свойств 
механических колебательных систем с двумя степенями свободы основаны на определении экстремальных 
значений частотных функций и функции демпфирования, определяемых из соотношений, связывающих 
кинетическую, потенциальную энергию и значения функции рассеяния энергии. Математические модели 
строятся на основе формализма Лагранжа, матричных методов, элементов теории функций комплексной 
переменной. 

Результаты исследования. Предложен метод построения частотных функций и функций демпфирования для 
класса механических колебательных систем с двумя степенями свободы на основе аналитических выражений, 
отражающих особенности соотношения потенциальной и кинетической энергии системы с учетом сил вязкого 
трения, представленных диссипативной функцией. Выведены общие аналитические выражения для частотной 
функции и функции демпфирования. Для механических колебательных систем с упруго-демпфирующими 
элементами, обладающими свойствами симметрии, проведен графоаналитический анализ экстремальных 
свойств соответствующих частотных функций и функций демпфирования. Представлены результаты 
численных экспериментов. Предложен критерий классификации частотных функций и функций 
демпфирования на основе топологических особенностей графиков соответствующих функций. 

Обсуждение и заключения. Разработанный метод построения частотных функций и функций демпфирования 
может быть использован для отображения динамических особенностей механических колебательных систем. 
Предложенный матричный метод построения частотно-демпирующей функции для системы двумя степенями 
свободы может быть распространен на механические колебательные системы, рассматриваемые в различных 
системах координат. 


Ключевые слова: механическая система, динамические связи, частотная функция, функция демпфирования, 
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Введение. Методам использования механических колебательных систем в качестве расчетных схем в 
задачах оценки динамических свойств технических объектов, работающих в условиях интенсивных 
вибрационных нагрузок, уделяется существенное внимание [1-9]. К числу распространенных подходов к 
оценке динамических свойств механических колебательных систем можно отнести методы, основанные на 
экстремальных свойствах отношения потенциальной и кинетической энергии [10, 11]. Методы, построенные на 
основе энергетических отношений, получили развитие в использовании частотной функции, как функции 
коэффициента форм связности координат механической системы, для оценки особенностей динамики 
механических колебательных систем без учета сил трения [12—15]. 

Вместе с тем, методы оценки динамических свойств механических колебательных систем, с учетом сил 
вязкого трения на основе частотной функции, требуют детализации представлений в зависимости от величины 
вязкого трения. Это связано, в частности, с тем, что для систем с апериодическим движением понятие частоты 
колебания может терять смысл. 

Предлагаемая работа посвящена развитию метода оценки свойств механических движений на основе 
развития понятия частотной функции, когда вводится дополнительная функция демпфирования, отражающая 
особенности учета сил вязкого трения в зависимости от коэффициентов форм связности. 

Материалы и методы. Рассматриваются свободные движения механической упруго-диссипативной 
системы с сосредоточенными параметрами с двумя степенями свободы. Принципиальная схема системы 
представлена на рис. 1. 


У! У? 


[2 Бо Ь? 


ы оС оС А 


Рис. 1. Механическая колебательная система с учетом вязкого трения 
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Обобщенные координаты у, у» обозначают смещение массо-инерционных элементов т, т, 


относительно положения статического равновесия. Кинетическая энергия ТГ, потенциальная энергия П и 
функция рассеяния Ё имеют вид: 


Г. | В 
Т=оту +5тьу, , (1) 
Ре И 2 
В. т м. > (2) 
1 В 1 : . 2 1 ‚2 
г И +5%0,-Я) о | (3) 
Система уравнений Лагранжа второго рода имеет вид: 
Ч ОТ ОТ ОоЕ 
+ =0; 
и, д, 9 7 
Я ОТ ОШ ОЕ 
ати 
@ 0, Оу, 0, 
После подстановки выражений Т,П, ЕЁ система дифференциальных уравнений (4) принимает вид: 
ы +(5 +5), Б.у, + (№ +А) у, — Ко, = 0; (5) 
т, у, + (6 +5, )у, -Б У + (& + )у, —К у, =0. 


Формы свободных движений представленной системы (5) в общем случае определяются собственными 
числами характеристического уравнения с учетом их кратности. Рассматривается случай простых корней. 
Таким образом, пусть решение у = у (Ру =у,() системы (5) представимо в виде: 


р”, (6) 


а - |И . . 
где у - — вектор-решение, У = —щ числовой вектор, р=о+ 76 — комплексный параметр, Е — 
2 2 


переменная времени. Полагается, что начальные условия согласованы с видом искомого решения (6). 

Задача заключается в построении и оценке экстремальных свойств функций, отображающих 
характеристики собственных движений системы с учетом сил вязкого трения. 

Результаты исследований 


1. Построение частотной функции и функции диссипации на основе энергетического 
соотношения. Система (5) в обозначениях (6) имеет вид: 


бы -№Р-№ | (1) 
-6Р-№ т, р” + (Б, +Ь, р + К, || № 
Введем обозначения: 
= а —Б, с" —№ | (8 
О т, —/ ВБ -Ь, —«ю Ю№-К, 
С учетом (8) матричное соотношение (5) принимает вид: 
(р?А+рВ+С)У =0. (9) 


На основе матричного соотношения (9) могут быть получены различные скалярные уравнения, с 
учетом экстремальных свойств которых, в свою очередь, могут быть определены свойства решения р уравнения 


(9). В частности, скалярное умножение левой и правой части равенства (9) на вектор У приводит к скалярному 


выражению: 
р? (А7,У)+р(В7,У)+(СТ,У)=0. (10) 
С помощью подстановки р=о+ 0 приведем (10) к виду: 
(0*—@7 +2/06)(А7,7)+(в+ /)(ВУ,У)+(СУ,У)=0, (11) 
Пусть для координат вектора У выполнено соотношение: 
У, =а7, (12) 
где © — коэффициент формы связности. В этом случае вектор У может быть представлен в виде: 
та. (13) 
где @ = В — вектор, определяемый коэффициентом связности @. После подстановки (13) выражение (11) 


принимает вид: 
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(9—0 +2/о)А, + (0+ /)В, +С, =0, (14) 
где А = (Аа, 6) ‚ В = (Ва, 6) ‚В. = (Ва, 6) — скалярные функции коэффициента формы связности &@. 
Уравнение (14) может быть представлено в эквивалентном виде: 

<’ А, =о”А, +оВ, +С.; 
20%А, +®В, =0. 


Найдём решение системы (15) относительно ©,©, как функций ©. Функция (0) рассматривается как 


(15) 


частотная функция, о(9) — функция демпфирования. Особенности решения системы определяются знаком 
2 
выражения о’4А +оВБ, +С.. 
Учет диссипации. Уровень диссипации энергии может быть охарактеризован дискриминантом 


В’ -4А,С,. При условии: 














В” <4АС, (16) 
понимаемом как малость трения, решение (15) может быть представлено в виде: 
2 С. В. 
А, 24)’ 
{0 {0 (17) 
В. 
о=- 
24. 
Необходимо отметить, что при выполнении условий малости сил вязкого трения (16) выполнено равенство: 
С 
2 2 
ю +о =—“. 18 
7 (18) 
При условии «больших сил вязкого трения»: 
В’ >4А.С., (19) 
решение (15) может быть представлено в виде: 
«=0, (20) в (@) 5. (21) 
— : о = — — 
| 2А, 
Функция демпфирования имеет две компоненты с’ (0), ©,(@,) такие, что: 
В С 
и, > (22) 


Таким образом, в зависимости от уровня сил вязкого трения, возможны различные виды представления 


частотной функции и функции демпфирования. Если В ^-4А.С, <0, то определена частотная функция > (0) 
и одна компонента функции демпфирования ©(0). Если В^-4А.С. >0, то предполагается, что частотная 
функция ©’(9) принимает нулевые значения, а функция демпфирования имеет две различные отрицательные 
компоненты с’ (©), с, (0). Что касается условия В *-4А.С, =0, то оно требует отдельного анализа. Вместе с 


2 
тем, условие В ^-4ААС =0 может быть интерпретировано как граница двух различных режимов движения 
[02 [02 [02 


механической системы. 

Представленные аналитические выражения частотной функции, функции демпфирования и условий 
«малых» и «больших» сил вязкого трения могут быть детализированы при рассмотрении частных вариантов 
механических колебательных систем, полученных на основе системы с двумя степенями свободы. 

2. Частотная функция и функция демпфирования для механической системы с двумя степенями 
свободы. Рассматриваются варианты параметров для механической системы, представленной на рис. 1. 
Предполагается, что множество граничных параметров, разделяющих режимы движение при малых и больших 
силах вязкого трения, определяется уравнением: 


В.’ =4А.С., (23) 
где: 
А =т +т, о”, (24) 
В .=(Ь,+Ь, ог —-2а6, +В. +В, (25) 
С =(® +К,)а” =? ++. (26) 


Условия малости сил вязкого трения имеют форму: 
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В.’ <4АС.. (27) 


В этом случае частотная функция <’ (0) и функция демпфирования о(0): 


2 
(а) = (®-А,) о — 29% ++ [1 (+, — 296, +В, +В, 
т, и 2 т, +т, а” (28) 
в 
ое Вы 
2 т, + т. 


Условия больших сил вязкого трения имеют форму: 
В’ >4АС.. (29) 


При условиях (29) функции ©’ и о(а) имеют вид: 


_ 1 (6+6, )7 —-2а6 +В. +В _ 


2 т, +т, а” 


О (0) — 








— [16 -+Ь)а” —296, +6, +6, (+6)? -20% +6 +8. 0) 
2 т, +т, о т, +т, о 
я 
_2 
И ав 
2 т, + та 


1 (6 -Ь, а —296, +В, +В, 


й 
Е т, + т, о” | 


Представленные выражения отражают движения в форме экспоненциального убывания при отсутствии 





_ (+) — 20, + + 


2 
т, + т.о 





колебаний. 
3. Особенности частотных функций и функций демпфировании для симметричных механических 
колебательных систем. Рассматривается механическая колебательная система с упруго-демпфирующими 


элементами, на значения параметров которых наложены условия симметрии в виде ВБ =Б, =ВБ =Б, 


К =А, =А =К. Принципиальная схема представлена на рис. 2. 


У? У? 
Ь Ь Ь 
т тТ>? 
у: К 
ОО О: © 
\\\ \\\ 


Рис. 2. «Симметричная» механическая система 


Система дифференциальных уравнений (5) имеет вид: 


ы +2 — 9, +25 №, =0; ых 
ту, +2Бу, -Бу +2Ку, -Ку =0. 
Функции А., В., С, могут быть представлены выражениями: 

А =т +т, о”, (32) 

В .=2Ы о’ -@+1), (33) 

С =2К(а/ -@-+1). (34) 


На основе представленных компонент могут быть построены частотная функция, функция 
демпфирования и сформулированы условия малости сил вязкого трения. 
Учет сил вязкого трения. Условие малости сил трения могут быть представлены неравенством: 


Во С 
В | <—. (35) 





[0 
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После подстановки функций (32)—(34) условие малости сил трения (35) может быть представлено в 
виде: 
у, <М(а), (36) 
.. й [т + т, а.” .. .. 
где у, = — — обобщенный вязко-упругий параметр, М (а) =--————— — обобщенный массо-инерционный 
4К 2 а -@а+1 
коэффициент, зависящий от коэффициента формы @&. График функции М. для каждого фиксированного 7, 
определяет множество значений ©, при которых выполнено условие малости сил трения. 


В качестве примера на рис. 3 приведен график параметризующей функции М(а). Функция М(@) 


2 


т 
обладает глобальным минимумом М, и максимумом М,, при а —00 М(@)-— Е 


Ус 
(1) 








— 


— — = — — —— -_ — ——_—_—_— =—— 


-10 -5 0. 5 10 





{0 


Рис. 3. Параметризующая функция М (0): 1 — уровень глобального максимума М, , 3 — уровень глобального минимума 


т 
М, , 2 — уровень горизонтальной асимптоты — 


Изменение параметра ‘у, в интервалах (,М.) (М.Р), С,М,) (М9) определяет характерные 


интервалы коэффициента формы ©, в которых выполнены условия малости сил трения. 
При малых силах трения частотная функция и функция демпфирования имеют вид: 


— 2к(о” —9+1 _ Ба” —@-+1) и 


о (©) т, + та т +т.а? |’ 
1 2 1 2 
(37) 
(о? — а +1 
= 
т, + т, а, 


При больших силах трения, при которых реализуется апериодическое движение системы, частотная 
функция равна нулю, а функция демпфирования имеет две компоненты: 








Ф’ =0; 
(о? — а +1 (а? о -+1 | 2(о? —а-1 
ыы. (88) 
т, + т. а т, + т. т, + т. 
(а? — “+1 а, ЭЕы а 
(а-я, |5 2% ) 





й 2 
т, + т.“ т, + т, о, 


На основе аналитических представлений частотной функции и функции демпфирования могут быть 
определены характерные варианты и особенности экстремальных свойств соответствующих функций с учетом 
сил вязкого трения. 

Обсуждение и заключение. Интерес представляют характерные варианты частотных функций и 
функций демпфирования в зависимости от условий малости сил вязкого трения. Рассматриваются варианты 
значений величины 7,, определяющие характерные интервалы коэффициента формы &@, в которых 
выполняются условия малости сил вязкого трения. 
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1. Пусть 7, Е (0, М,). Рассматривается пример 7, = 0,1. В этом случае выполнены условия малости сил 
трения для любого коэффициента формы &е (-—0..оо). На рис. 4 и 5 представлены частотная функция 6’ (0) и 
функция демпфирования ©(0) для механической упруго-диссипативной системы с параметрами Б=1, т =3, 
т, =4,К=З. 





-0.5- в 
Рис. 4. Частотная функция ©’ (0): Рис. 5. Функция демпфирования ©(0,): 
(1) и (2) — экстремальные уровни (1) и (2) — экстремальные уровни 


Корни уравнения, представляющего собой равенство нулю соответствующего определителя 
Ар? +Вр+С|=0, (39) 
составляют р, =®,+10,, 1=1..4, где в = 0,91; ©, =1,56; ®, =-1,56; в, =-0,91; в =-0,14; в, =-0,44 ; 
с, = -—0,44; с, =-0,14. На рис. 4 частотная функция достигает экстремальных значений, равных квадратам 
частот ©,‘ 22,46 и в* = 0,82. На рис. 5 функция демпфирования достигает экстремальных значений, которые 
составляют ©; =-0,44 и оа,=-0,14. Свои экстремальные значения частотная функция и функция 


демпфирования достигают, когда коэффициенты формы составляют © * = -0,65 и ©,* =1,15. 


т 
2. Пусть у, Е(М,,—). Рассматривается механическая система с параметрами Б=4, 7, =1,33. На 

2 
рис.6 и 7 представлены соответствующие частотная функция и функция демпфирования. Множество 
коэффициентов форм, для которых выполняется условие малых сил трения, составляет: (—00, 0 ) (4,00), где ©, 
— корни уравнения М(о)=у,. Для параметров Р=4, т=3, т, =4, К=З3 корни характеристического 
уравнения (39) имеют вещественные части, представляющие собой коэффициенты диссипации, ©, = -—0,57; 
с, =—1,08; а, =-2,46; <, =-0,57 и мнимые части, представляющие собой частоты, ©, = 0,73; ©, =0; ®, =0; 


> —. 
©, = —0, 73 . На рис. 6 на интервале (—, © )\/(©,,с0) частотная функция имеет локальный минимум в,’ = 0,53. 


у. 
@ [9 


0.6 





(1) 
(2) 
-31 б 
9:1 
Рис. 6. Частотная функция : Рис. 7. Функция демпфирования : (3) -экстремальный 
(1) — экстремальный уровень в точке о. уровень в точке о ; (1), (2) — экстремальные уровни в точке 


* 
О на интервале двузначности 
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На интервале (©,о,) частотная функция равна нулю. В свою очередь, на рис. 7 функция 
демпфирования в интервале (0,0) двузначна и достигает одновременно двух экстремальных значений 
с, =—1,08 и о, =-2,46 в точке а,*. В области (—°,0 )\/(©,,00) функция демпфирования однозначна и 
имеет один локальный экстремум с, = -0,57 в точке @, *. 

3. Пусть уу Е М › ). Рассматривается механическая система с параметрами Б=6; т =3; т, =4; 

2 
К =3. Корни характеристического уравнения (39) имеют вещественные в, =-0,85; в, =-—0,56; в. я 4,74; 
с, =-0,85 и мнимые части ©, = 0,36; ®, =0; ®, =0; <, =-—0,36. Условия малости сил трения выполняются 
на интервале (0,0), где о = 0,63; @, = 2,37. На рис. 8 соответствующая частотная функция положительна 
только на интервале (9,0). Локальный экстремум частотной функции равен 6.” = 0,13. Вне интервала 
(0, %,) частотная функция равна нулю. 

Функция демпфирования на рис. 9 в интервале (—0,0/) двузначна и достигает одновременно двух 
экстремальных значений в точке а, *, которые составляют а, = 4,74 и в, =-0,56. В интервале (&,, оо) 
функция демпфирования также двузначна. В интервале (0,0) функция демпфирования однозначна и имеет 


один локальный экстремум о, = 0,85 в точке ©, *. 





(1) [9 
-4 -2 0 2 4 
[\, 
„1 —ы-{ 
(2) 
_-3 
| (3) 
] 
_-5 б 


в Рис. 9. Функция демпфирования: (1), (3) — экстремальные 


. * 
Рис. 8. Частотная функция: уровни в точке 0, на интервале двузначности, 


* 
вы к * 
(1) — экстремальный уровень в точке 0, (2) — экстремальный уровень в точке (©, 


4. Пусть у, Е (М, ). Для параметров Б = 6,52; т =3; т, =4; К=3 характеристическое уравнение 
(39) имеет только вещественные корни © =-0,50; в, =-0,88; в, =-0,96; в, =-5,26. Соответствующая 


частотная функция и функция демпфирования представлены на рис. 10 и 11. Интервал выполнения условий 
малости сил трения вырождается в пустое множество. 





(2) (1) (1 
> 
у 
в. 
(3) 
=Э 
(4) 

-4 
3 
-5 = 
ы 
(>) 
-6] св х 
Рис. 10. Частотная функция: Рис. 11. Функция демпфирования, образованная двумя ветвями: и 2 

вырожденный случай (1), (2) — экстремальные уровни «верхней» ветви в точках 0 и ©. ; 


* * 
(3), (4) — экстремальные уровни «нижней» ветви в точках 0 и © 
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На рис. 10 частотная функция равна нулю на всей числовой оси. На рис. 11 соответствующая функция 


демпфирования двузначна на всей числовой оси и имеет 4 локальных экстремума о, =-0,50; о, =-0,88; 


нЕ НЕЕ 


= 


* * 
с, 0,96; <, =-—5,26 в точках @ и а, 


(д) нь 








Рис. 12. Пиктограммы топологических особенностей графиков частотной функции и функции демпфирования: 
(а)-(г) — пиктограммы, сопоставленные графикам частотных функций, представленных на рис. 4, 6, 8, 10 соответственно; 
(9)-(3)— пиктограммы, сопоставленные функциям демпфирования на рис. 5, 7, 9, 1] соответственно 


Представленные на рис. 12 пиктограммы, сопоставленные графикам функций на рис. 4—11, отражают 
ряд топологических особенностей графиков частотной функции и функции демпфирования. К таким 
особенностям можно отнести форму графика в виде одной кривой, наличие на графике точек раздваивания 
одной кривой на две, наличие двух непересекающихся кривых или «кольца». 

Таким образом, достигаемые экстремальные значения частотной функции и функции демпфирования 
связаны с динамическими характеристиками механической колебательной системы с учетом сил трения. В 
частности, экстремальные значения построенных частотной функции и функции демпфирования связаны с 
собственными частотами и диссипативными коэффициентами затухающих колебаний. При этом вопрос 
существования экстремальных значений частотной функции и функции демпфирования, не совпадающих с 
квадратами частот собственных колебаний, требует дополнительного рассмотрения. Вместе с тем, можно 
предположить, что формы частотных функций и функций демпфирования, обуславливающие формы 
свободных движений механических колебательных систем с трением, представляют интерес для оценки более 
широкого круга динамических свойств. 

В плане практической реализации возможного управления колебательными режимами механических 
систем на основе коэффициента связности не существует принципиальных препятствий. К примеру, возможные 
системы управления динамическим состоянием могут включать датчики амплитуд колебания в контрольных 
точках вибрационной технологической машины или транспортного средства. Вместе с тем, вопросы 
построения подобных систем требуют детального рассмотрения широкого круга особенностей, связанных с 
техническим объектом. 

В заключение в качестве выводов представленных исследований можно отметить следующие моменты. 

1. Для механической колебательной системы с учетом сил вязкого трения разработан метод построения 
частотной функции и функции демпфирования, зависящих от коэффициента форм связности координат 
свободных движений. Показано, что совокупность экстремальных значений частотной функции и функции 
демпфирования отображает собственные характеристики упруго-диссипативной механической колебательной 
системы. 

2. Показано, что частотная функция и функции демпфирования для механической колебательной с двумя 
степенями системы с учетом вязкого трения может быть представлена в двух вариантах, определяемых условиями на 
величину сил вязкого трения для фиксированного коэффициента форм связности; для условий малых сил вязкого 
трения значения частотной функции принимают положительные значения, а функция демпфирования имеет одну 
отрицательную компоненту; если выполнены условия больших сил вязкого трения, то частотная функция принимает 
нулевые значения, а функция демпфирования имеет две отрицательные компоненты. 

4. Предложен метод построения возможных вариантов частотных функций и функций демпфирования для 
различных значений параметров системы на основе параметризующей функции, позволяющей определить области 
значений коэффициента форм связности, в которых выполняется условие малости сил вязкого трения. Предложен 
критерий классификации частотных функций и функций демпфирования в зависимости от топологических 
особенностей их графиков. 

5. Матричный метод построения частотно-демпирующей функции для системы двумя степенями 
свободы может быть распространен на механические колебательные системы, рассматриваемые в различных 
системах координат. 
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6. В качестве физической интерпретации коэффициента связности, используемого в частотной функции 
и функции демпфирования, может рассматриваться рычажная связь в форме передаточного отношения, 
выраженного через отношения амплитуд колебаний координат парциальных блоков. Рассматриваемое 
отношение, наравне со статическим состоянием, может быть определено для режимов установившихся и 
затухающих колебаний. Таким образом, развивается концепция, в рамках которой отправной точкой анализа 
механической системы является рычажная связь. 
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